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In this work we calculate the meson loop effects on the pion electromagnetic form factor and 
perform a comparison with experimental data. We show that, even though the meson cloud is not 
the only nonperturbative process to be considered, its contribution is significant. 
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I. INTRODUCTION 

The pion electromagnetic form factor has been under intense investigation both experimentally and theoretically. 
Yet there are still some open questions. Perhaps the most intriguing one is: why experimental data in the region of 
Q 2 > 2 GeV 2 are still so far away from the perturbative QCD (pQCD) expectation? Or equivalently, what is the 
nonpertubative dynamics responsible for this behavior? Since long ago it has been suspected that the transition to 
the pQCD regime might happen at intermediate values of Q 2 (5 < Q 2 < 15 GeV 2 ) [1-3]. However, the data from 
Jefferson Laboratory (JLAB) [4] are not only far from the pQCD prediction but follow a different trend. 

Lattice calculations, in spite of their rapid progress [5], are not yet in position to answer the questions above. 
In order to elucidate the subject many models have been advanced. A variety of constituent quark models has 
| been considered [6-9]. Most of them are based on solutions of the light front Bethe-Salpeter equation with some 
approximations . In general, they achieve good agreement with data. Along a different line, the instanton model 
Q\ ' developed in Ref. [10] was able to give a reasonable description of data in the region Q 2 > lGeV 2 . In a similar 
CO . approach the same conclusion was found in [11]. 

In this work we consider the contribution of the meson cloud around the pion to its eletromagnetic form factor. 
This nonperturbative, hadronic effect was shown to play an important role in the electromagnetic form factor of the 
nucleon [12-16]. Recently the effect of the pion cloud was considered in the study of the pion form factor in the 
timelike region [17]. The p meson, provenient from the photon was converted into a virtual pion loop. Using the 
same approach and the same kind of vertices, we compute mesonic loop contribution to the pion form factor in the 
spacelike region. 

The idea of describing a pion with cloud states has already been used in [18] and also by us previously in [19], where 
' we have studied D~ / D + asymmetry produced in irp — > DX collisions. In that work we could see that the inclusion 
^ , of the meson cloud leads to a good agreement with experimental data. 

In the next section we will review a few formulas used in the form factor calculations and define the dominant 
processes. In section III we discuss the inputs for the calculations and the numerical results. In section IV we present 
some concluding remarks. 
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II. LOOPS AND THE PION FORM FACTOR 



The pion electromagnetic form factor is defined by: 

< n-(p')\J u \Tr-(p) >=e(p + p') t ,F w (q 2 ), q = p'-p, (1) 

where is the electromagnetic current. As in the case of the nucleon, we can represent the pion as a "bare pion" 
1 7r >o which can fluctuate into virtual Fock states composed by a pseudo-scalar (P) meson and a vector (V) meson, 
such as: 

|7T > = Itt" > +|7T-p° > +|ttV > +\K-K°* > +\K°K~* > ... (2) 

Thus, we can understand the photon-pion interaction as being an interaction between the photon and the constituents 
of the cloud. Usually the meson loop contributions are supposed to be dominant at low values of Q 2 , where the 
incoming photon does not yet resolve quarks inside the pion. At higher values of Q 2 , quark pointlike structures 
should become visible and a description of the problem in terms of hadronic degrees of freedom should loose validity. 

However, due to the present difficulty in understanding the intermediate Q 2 data, we consider the meson cloud 
contributions up to Q 2 = 6 GeV 2 . 
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In Fig. 1 we show the first terms generated by the expansion (2): first the interaction between the photon and 
the bare pion (lo) and then the interaction between the photon and states composed by a pseudoscalar meson and a 
vector meson \PV > (la-d). The full and open circles represent the dressed form factor and the form factor of the 
bare pion, respectively. 

Before presenting the calculations some remarks are in order. Quite generally, the structure, or size, of any hadron 
depends on its quark content, on its meson cloud and on its off-shellness. We want to isolate and compute the strength 
of the meson cloud contribution. Looking at Fig. 1 and comparing the pion in the l.h.s. of the equality with the pion 
in a) we realize that the first one is on-shell, whereas the second one is off-shell. Therefore they are not exactly the 
same object and do not have the same structure. Moreover the pion in Fig. la) is a bare pion whereas the pion on the 
l.h.s. of the equality is composed by a bare pion plus the sum over all cloud corrections. With these considerations 
we want to emphasize that the l.h.s. pion and the pion within the loop are not the same. 
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FIG. 1. Photon-pion interaction in the meson cloud approach. 



In order to calculate de diagrams of Fig. 1 we start from the interaction Lagrangian given by: 



(3) 



where P and V denote a pseudoscalar and a vector meson respectively. From £ we derive the Feynman rules with 
which we write the vertex function associated with diagram (a) of Fig. 1 : 



ri a) (p,p') = ig 2 Q P J ^F(k)F(k + q)A(k)A(k + q)p' D^(p-k)p a (2k + q) fl 
and with diagram (b): 



(4) 



(p + p' - 2fc) p gg v - (p' - k) v g 6fl - (p - k) e g nil 



D a \p-k) Pa 



(5) 



where Qp and Qy are the charges of the pseudoscalar and vector mesons, and A(fc) and D al3 (p - k) are their 
propagators. They are given by: 



D af3 (p-k) = 



g afi - {p-k) a {p-kf /m v 

(P — k) 2 — ffiy 



and A(fc) 



1 



k 2 
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where mp and my are the masses of the pseudoscalar and vector mesons respectively. 

In order to account in some way for the finite extent of the mesons appearing in the loops of Fig. 1, we have 
included form factors in the hadronic vertices. For simplicity we have chosen a monopole form: 

m = («) 

where A is the cut-off parameter. Although there is no rigorous justification for this choice, form factors of this type 
have been used for example in Ref. [17]. These form factors render all the following loop integrals finite. 



In the presence of the electromagnetic field the nonlocal mesonic interaction of Eqs. (3) and (6) gives rise to 
vertex currents. In order to maintain gauge invariance we introduce the photon field by minimal substitution of the 
momentum variable in the form factors. This procedure generates nonlocal seagull vertices (s) [20]. They are: 



for diagram (c) and 



= 9Qpp' a (<! 2% (fc _^ ( 2 fc) _ A2 (7) 



*U = gQpp a (q + 2k), {k + F ^ ] _ A2 (8) 



for diagram (d). Using the expressions above to write the T,(p,p') functions for diagrams (c) and (d) we find 
respectively: 



Tf°{p,p')=ig 2 Q P / ^ A(fc) F 2 (fc) (fc _ -1 _ A2 p' (j £>°V - k) P ' a (q 2k), 



(9) 



and 



r^(p,p') = - ig 2 Q P J -0^ A(k) F\k) - - -1 _ A2 Pp D"P(p - k) Pa (q + 2k), (10) 

Due to the presence of the derivative in Eq. (3), the minimal substitution also generates current-vertex couplings 
(v) even if the mesonic form factors are not present. They are: 



v c = 9 ( Qd F(k -q) + Q v F(k) ) (11 ) 

for diagram (c) and 

v d = -9 [ Qd F(k + q) + Q v F(k) ) (12) 



for diagram (d). Repeating the procedure used above, we write the T,(p,p') functions for diagrams (c) and (d) using 
(11) and (12): 

r^(p,p')=ig 2 J A W F{k)PpD^{p' -k)^Q P F(k-q) + Q v F(k)^ (13) 

and 

T^(p,p') = zg 2 J ^A(fc) F(k)p D^(p-k)^Q P F(k + q) + Q v F(k) j (14) 
In order to obtain the total vertex function T„ we add (4), (5), (9), (10), (13) and (14) to find: 

v,{p, P ') = F l {Q 2 )p ll + F 2 {Q 2 )p' fl 

The pion electromagnetic form factor is associated with the term: 
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F+(Q 2 ) = 



F^Q 2 ) + F 2 {Q 2 ) 



(15) 



Since in the case of the pion it has not been discussed before we shall, in the remainder of this section, briefly show 
that the inclusion of the terms (9), (10), (13) and (14) is crucial for the fullfilment of the Ward-Takahashi (W-T) 
identity. The W-T identity reads: 



where 



^ = Q^s(p) - E( P ')) 

H{p) = -ig 2 J ^ A{k)F 2 (k)p fi D^{p-k)p a 

F(k)F(k + q)p' D^(p-k) Po 



A(k)-A(k + q) 



For diagram (a) we have: 

a^rW = -ia 2 Op ( d * k 
and for diagram (b): 

q^ = -ig 2 Q v J ^ A(k)F 2 (k)p' p 
For the seagull coupling of diagrams (c) and (d) : 

q^=-ig 2 Q P J -0^ A(k)F(k)x 

-p'pD^ip' -k)p' a (F(k) - F(k - q)j + Pp D a f } (p-k) Pa (F(k)-F{k + q)j 



D af3 (p' -k)-D af \p-k) 



(16) 



(17) 



(18) 



(19) 



and for the vertex coupling of diagrams (c) and (d): 

d 4 k 
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(2tt)4 
A(k)F 2 (k) 



A(k)F(k) 



F(k - q)p' fj D ati (p' -k)q a - F(k + q)p a D afj (p - k) qp 



p' & D a(3 (p' -k)q a - p a D af3 (p - k) qp 



(20) 



Adding the contributions (17), (18), (20) and (19), it can easily be seen that: 

d 4 k 



q»T ll = -ig 2 Qp 



-ig Q v 



(2tt)* 
d 4 k 



A{k)F 2 (k) 
A(k)F 2 (k) 



(2tt)< 

= Q P (S(p) - E(p')) + Qv (s(p) - E(p') 
Due to charge conservation, Qp + Qv = Qtt, and thus we have: 

g M r^ = Q w ^s(p)-E( P ')) 



PpD^(p - k) Pa -p'pD^ip* - k)p' a 
P fi D a P{p -k)p a - p'pD^{p' - k) p' Q 



(21) 



As it can be seen from the last factor in expression (4), we are using a pointlike coupling in the virtual pion - photon 
- virtual pion (Fig. la). Since this pion has already some structure, we might introduce a form factor in this coupling, 
given by its quark structure (already present in the bare pion). However, we are interested only in the meson cloud 
and not in modelling the bare pion. The most simple and consistent procedure is to take, as we did, a pointlikc 
coupling and keep in mind that this will maximize the contribution of the cloud, since any possible form factor would 
be a decreasing function of Q 2 , with maximum value equal to one (at Q 2 = 0). In this way we are computing the 
upper value for the cloud contribution. 

Being restricted to the cloud contribution and making no attempt to accurately describe the data, we have only 
one parameter (the cut-off Lambda in our meson-meson-meson form factors), which is constrained by other types of 
phcnomcnological analyses. We have thus little freedom and therefore some predictive power. 
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III. NUMERICAL RESULTS AND DISCUSSION 



In order to compute the form factor (15), it is necessary to know the value of the relevant masses, coupling constants 
(g) and cut-off parameters (A). The g T7rp , g^KK* coupling constants have been measured experimentally [21]. As the 
masses are all known, only the cut-off values remain as free parameters. So far there is no systematic study which 
indicates the most appropriate cut-off value for the vertices involving pion-meson-meson fluctuations. Since we are 
using a monopole parametrization for these vertices, we have chosen cut-off values close to the corresponding vector 
meson masses. In Table I we present all the values used in these calculations. 
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Table I: Values of coupling constants and cut-off used in the vertex functions. 

In order to explicitly show the relative strength of each cloud state, we show in Fig. 2 an example of result obtained 
with (15) and a cut-off fixed to A = 0.9 GeV GeV for the sum of \np > states (dashed line), the sum of \KK* > states 
(dot-dashed line), and the sum of all of them (solid line). We include the experimental data from the F„ collaboration 
[4], from Amendolia et al. [22], from DESY [23] and from Bebeck et al. [24]. As it can be seen, all the curves rise , 
reach a maximum and then drop continuously. The heigth and relative strength of each contribution is controlled by 
the masses, couplings and cut-off values. 
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FIG. 2. Contribution of the meson cloud states to the pion form factor. The black circles, crosses, empty circles and black 
triangles are the experimental data from Refs. [4], [22], [23] and [24] respectively. In this figure we used A mTp =A J rKK * =0.9 GeV. 
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FIG. 3. Contribution of the sum of all meson cloud states to the quantity Q 2 F W (Q 2 ), with different cut-off choices. In this 
figure A = A^p = AkKK * ■ 



In Fig. 3 we show the contribution of the sum of all the meson cloud states to the pion form factor for several cut-off 
choices. We can observe a strong cut-off dependence. Nevertheless it is possible to conclude that the cloud contribution 
is significant specially in the region 1 < Q 2 < 6 GeV 2 . Compared to the case of the nucleon electromagnetic form 
factor [13] the meson cloud plays here a less important role. 

For each choice of cut-off parameter in Fig. 3 we have calculated the corresponding mean square radius < r 2 >, 
according to the usual definition: 



< K >= 



-6 



dF^(Q 2 = 0) 
dQ 2 



(22) 



The results are presented in Table II. 
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Table II: The mean square radius < r 2 > for different cut-off values. 

We can see that, whatever the cut-off choice is, the mean square radius is much smaller than the experimental value 
< rl >= 0.439 ± 0.008 fm 2 [22]. 

In Fig. 4 we plot the form factor F+(Q 2 ) without multiplying it by Q 2 to show in detail the behavior of our curves 
at low values of Q 2 . 
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FIG. 4. The same as figure 3 for F W {Q 2 ). 



The results in Figs. 3 and 4 indicate that the cloud contribution underestimates the data in the low Q 2 region. At 
higher values of Q 2 (Q 2 > 2 GeV 2 ), it starts to be closer to data. In the asymptotic Q 2 limit the cloud contribution 
vanishes. Since we did not make any model for the bare pion, its influence on the form factor is a priori unknown. 
Looking at the results presented in Fig. 4, we can just observe that there is room for a significant bare pion 
contribution. 

Perturbative QCD predicts that in the high Q 2 region the form factor should be much smaller than the actually 
measured data. The cloud component is then a candidate to be the nonperturbative complement to pQCD. In order 
to further investigate the relation bewteen the cloud and pQCD, we make now the simplifying assumption that in 
the asymptotic Q 2 region the bare pion contribution is simply given by pQCD. The physical picture becomes quite 
simple: a highly virtual photon probes a pion which is either the "core" or the cloud. In the first case the photon 
interacts with the quarks in the perturbative regime. In the second case, the photon interacts with a hadron of the 
cloud. In practice, the above assumption means adding to the cloud a contribution which is known and given by [10]: 

FAQ") = (23) 

where / ff = 92.4 MeV and a s = 0.4. This expression is plotted in Fig. 5 with a dotted line and it is added to the 
previously computed cloud contribution for A = 0.7 GeV, represented by the dashed line. We can see that the sum 
of the two contributions (solid line) gives a good description of data in the region Q 2 > 1 GeV 2 . In Fig. 5 we can 
clearly observe how, at large Q 2 , the cloud contribution dies out while pQCD takes over. Although our results are 
not definitive, they strongly suggest that a hard cut-off (A > 0.7 — 0.8 GeV) is hard to accomodate, if one is bound 
to reproduce data leaving room for a sizeable bare pion (pQCD) component. 
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FIG. 5. Fk(Q 2 ) with cloud (A = 0.7 GeV) and pQCD component. 



IV. CONCLUDING REMARKS 

Meson loop effects are a necessary consequence of quantum field theory and must be considered in the analysis of 
any hadronic observable. We have performed a straightforward extension of standard one loop calculations to the 
pion form factor. The main goal was to estimate the magnitude of the meson cloud contribution to F 7T (Q 2 ). It turned 
out to be less important than the meson cloud contribution to the nucleon form factor [13]. This is not unexpected 
since, because of the masses and corresponding virtualities, it is much easier for a proton to be in a neutron-pion state 
than it is for a pion to fluctuate into a pion-rho state. 

The size of the cloud contribution must also be cross checked with that in deep inelastic scattering (DIS), as discussed 
in [18]. In Fig. 4 we can see that the size of the cloud contribution depends on the value of Q 2 . At Q 2 around zero it 
represents something between 0.15 and 0.40 of the total form factor. This fraction is consistent with other estimates 
[18], coming from DIS and high energy inclusive hadron-hadron collisions. In [18], a smaller percentage was found for 
the cloud contribution. This discrepancy is probably due to the different methods used to extract the strength of the 
cloud states and also because different experimental data were considered. As Q 2 increases, the relative contribution 
of the cloud increases, but this region it is difficult to draw firm conclusions because the model gradually runs out of 
its domain of applicability and also, as stated above, our procedure gives an upper bound for the cloud contribution. 

The calculation presented here is model dependent in the sense that we have truncated a series (of cloud states) 
and replaced the effects of the higher order terms by form factors, which also account for the quark substructure 
of mesons. Since we have used the same type of form factor and the same cut-off for all vertices we have only one 
parameter, A, which was already constrained by our previous knowledge of hadronic interactions. It should lie in the 
interval 0.5 < A < 1.0 GeV. Our results confirm this expectation and also suggest that the cut-off should rather be 
in the range 0.6 — 0.7 GeV if we believe that that the cloud contribution is of the order of 15% or less. 

The fact that meson loops can not describe data in the low Q 2 region, regardless of cut-off choice, confirms that 
another dynamics is missing in this region, as, for example, the interaction of the photon with constituent quarks. 
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